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Abstract 

We consider a model of queues in discrete time, with batch services and arrivals. The 
case where arrival and service batches both have Bernoulli distributions corresponds to 
a discrete-time M/M/l queue, and the case where both have geometric distributions has 
also been previously studied. We describe a common extension to a more general class 
where the batches are the product of a Bernoulli and a geometric, and use reversibility 
arguments to prove versions of Burke's theorem for these models. Extensions to models 
with continuous time or continuous workload are also described. As an application, we 
show how these results can be combined with methods of Seppalainen and O'Connell to 
provide exact solutions for a new class of first-passage percolation problems. 

1 Introduction 

We consider a model of queues in discrete time, with batch services and arrivals. At the 
beginning of time slot n, there are X n customers in the queue. A number A n of customers 
then arrive, increasing the queue length to X n + A n . After this an amount S n of service is 
available, so that D n = mm(S n , X n + A n ) customers depart from the queue. Typically, we 
assume that the sequences A n and S n are random. 

This model has been studied in various contexts (sometimes described as a storage model 
rather than a queue |3J). There is a close correspondence between this model and another 
type of queueing model, in which the data S n represent inter-arrival times between successive 
customers and the data A n represent service requirements of customers. 

The case where the sequences S n and A n are independent and both consist of i.i.d. 
Bernoulli random variables corresponds to an M/M/l queue in discrete time. The case 
where the Bernoulli distributions are replaced by geometric distributions has also been stud- 
ied previously HHHU]. 

We generalise these two situations to the case where the distributions are a product of a 
Bernoulli and a geometric. 

We show that for appropriate choices of the parameters, the queue is reversible in equi- 
librium, and that various forms of Burke's theorem hold: for example, the departure process 
has the same distribution as the arrival process, and the queue-length at a given time is 
independent of the departure process before that time. The stationary distribution of the 
queue-length is also given by a product of a Bernoulli and a geometric. 

These properties make it easy to describe the stationary behaviour of several such queues 
in tandem. As an application, we show how this can be used to calculate rates of growth for 



certain first-passage percolation problems. This uses techniques developed by Seppalainen 
[T2"] and O'Connell [TO], and extends the class of "exactly solvable" first-passage percolation 
models. 

In Section [2] we describe the basic queueing model in more detail. In Section [3] we collect 
some definitions regarding probability distributions, in particular the distribution obtained 
by multiplying a Bernoulli and a geometric. 

In Section [4] we state and prove the reversibility results and versions of Burke's theorem, 
for queues with Bernoulli-geometric arrival and service processes. We discuss related results 
concerning systems of queues in tandem. We also note results about the stationary distribu- 
tion of the queue-length in the case where only the arrival process is assumed to have this 
form, and the service process consists of any i.i.d. sequence. 

As well as extending the Bernoulli and geometric cases, the Bernoulli-geometric model 
has another useful feature, that by taking a limit as the parameter of the Bernoulli tends 
to one can easily arrive at various continuous-time models in which arriving customers 
or offered service occur at times corresponding to points of Poisson processes. Taking an 
alternative limit, one can also consider models with continuous workload, where arrival and 
service batches are exponential rather than geometric. These extensions are indicated in 
Section PI 

In Section[5]we describe the application to first-passage percolation models, and give exact 
expressions for some time-constants. 



2 Queueing model with batch services and arrivals 

We describe the main queueing model of the paper. 

The queue is driven by an arrival process (A n , n £ Z) and a service process (S n , n € Z). 

At time-slot n £ Z, A n customers arrive at the queue. Then service is available for S n 
customers; if the queue-length is at least S n , then S n customers are served, while if the queue 
length is less than S n then all the customers are served. 

Let X n be the queue length after the service <SVi-i, before the arrival A n . From the 
description of the queue above, we have the basic recurrences 

X n +1 = [X n + An — S n ] + (2-1) 

(where denotes max{a;,0}). Similarly if Y n is the queue length after the arrival A n and 
before the service S n , then Y n — X n + A n , and 

Yn+l = [Y n — S n ] + + A n+ i. (2-2) 

In this paper we will almost always consider the case where the arrival and service processes 
are independent, and (i„,n£ Z) and (S n ,n € Z) are both i.i.d. sequences, with EA n < 
ES„. In this case (in fact, much more generally) we can define the queue- length sequence 
(X n ,neZ) by 

n-l 

X n = max V (A r - S r ) (2.3) 
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Figure 2.1: The evolution of the queue with batch services and arrivals 



(where a sum from n to n — 1 is understood to be 0) . This quantity is almost surely finite 
since the common mean of the S n is larger than that of the A n . Then the sequence (X n ) 



satisfies the recurrences (2.1), and that (X n ) is a stationary Markov chain. 

Let D n be the number departing from the queue at the time of the service S n . So 

D n = min(y n ,5 n ) 
= Y n — X n+ \ 
= X n + A n — X n+ i 
= Y n + A n+ i — Y n jf-\. 

Let U n be the unused service at the time of the service S n ■ So 

U n — S n D n 

= [S n — Y n ] + . 

See Figure [23] for a representation of the evolution of the queue along with its inputs and 
outputs. 

We also introduce some further quantities whose interpretation is ostensibly less natural 
(but see Section 2.1 1. Write I n = U n + A n+ i for the unused service plus next arrival, and 
T n = U n + A n for the unused service plus previous arrival. 

Note that although we have talked in terms of "numbers of customers" , there is no reason 
why the variables have to take integer values. We will also consider the case where A n and 
S n are non-negative real- valued random variables; here one might talk of "amount of work" , 
say, rather than "number of customers" . 

This model of a discrete-time queue with batch services and arrivals has been considered in 
various contexts, for example by Bedekar and Azizoglu [2], Ganesh, O'Conncll and Prabhakar 
4 , Draief, O'Conncll and Mairesse [3], and O'Connell [10]. 

Models that have been studied include the case where arrival batches and service batches 
are Bernoulli distributed, and the case where both have a geometric distribution; see the 
next sections for further details. The main result in this paper is to extend the versions of 
Burke's theorem obtained in these cases to the case of a distribution which is the product of 
a Bernoulli and a geometric. 
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2.1 Dual queueing model 

Our main batch queueing model is closely related to an alternative (and in fact more widely- 
studied) model of a single-server queue with first-in-first-out service discipline. 

This dual model uses the same variables and recurrences as above, but with different 
interpretations. Now A n represents the amount of time required by the nth customer for 
service, and S n is the interarrival time between customers n and n + 1 . 

Let X n be the waiting time of customer n between his arrival at the queue and the start 



of his service. Then X n obeys the same recurrence as at (2.1 1 above. Now Y n is the total time 
spent by customer n in the queue (including service), and D n is the time spent by customer 
n at the back of the queue. U n is the idle time of the server between departure of customer n 
and arrival of customer n + 1 , /„ is the interdeparture time between customers n and n + 1 , 
and T n is the time between the starts of service of customers n and n+ 1. 

See Draief, Mairesse and O'Connell [3] for extensive discussions of the relations between 
the two models. 

2.2 Geometric case 

We assume that the processes (A n ,n G Z) and (S n ,n G Z) are independent, and that each 
process is an i.i.d. sequence. 

Suppose now that both A n and S n have geometric distribution, with ¥(A n = k) = a(l — 

a) fe_1 and P(S n = k) = /3(l — /3) fe_1 for k = 1,2, For stability of the process, we require 

f3 < a. 



In this case, the single-server queue of Section 2.1 is an M/M/l queue in discrete time. A 
version of Burke's theorem was proved for this model by Hsu and Burke in 0. Among other 
properties one has that the arrival and departure processes have the same law, which (in the 
notation of Section 



2.1| says that (S n , «eZ) = {I n ,n G Z). 



On the other hand, Bedekar and Azizoglu [2] showed an analogous input-output theorem 
for the batch-queueing model; namely, that (A n ,n G Z) = (D n ,n G Z). 

These results arc unified by Draief, Mairesse and O'Connell 3 , who obtain a, joint Burke's 
theorem for the two models, namely that 

((A n ,S n ),n G Z) i ((£>„,/„), n G Z). (2.4) 

The same result also applies if the distributions of A n and S n are exponential rather than 
geometric. 

2.3 Bernoulli case 

Suppose instead that A n and S n both have Bernoulli distribution, with p — ¥(A n = 1) = 
1 - ¥(A n = 0) and q = ¥(S n = 1) = 1 - P(5„ = 0). For stability of the process, we require 
P < q. 

In our main batch-queue model, this means that each arrival batch is either empty or 
contains a single customer, and at each slot the available service is either 1 or 0. 

In this case, the batch-queue model in fact corresponds to an M/M/l queue in discrete 



time (hence to the dual queueing model in the geometric case of Section 2.2 1, since the 
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intervals between arrivals of successive customers are geometric (with mean 1/p) and the 
service time of a customer once he reaches the front of the queue is also ge omet ric (with mean 

corresponds to 



2.2 



l/q). So the statement that (S n ) — (I n ) in the geometric case of Section 
the statement that (A n ) = (D n ) in the Bernoulli case. 

Again this can be extended to a sort of joint input-output theorem for the two models in 
the Bernoulli case; namely one has that 



((A n ,S n ),n e Z) = ((D n ,T n ),ne Z). 
This result is contained in the proof of Theorem 4.1 of Konig, O'Connell and Roch [5]. 



(2.5) 



Note the difference between (2.4 1 and (2.5). In the geometric case, the output theorem for 



the single-server queue involves the interdeparture process (/„), while in the Bernoulli case it 
involves instead the process (T n ) of intervals between starts of service. One can see that I n 
is not Bernoulli in the Bernoulli case (indeed, one may have /„ = 2); also, one can see that 



T n and T„+i are not independent in the geometric case so that (2.5 1 certainly does not hold 
there. 



In Theorem 4. 1 we'll show that the common part of ( 2.4 1 and ( 2.5 ) , namely the result that 



(A n ) = (D n ), extends to a class of cases where the distributions of A n and S n are products 
of a geometric distribution and a Bernoulli distribution. 



3 Bernoulli-geometric distribution 

We first introduce some notation. X is said to have Bernoulli distribution with parameter p 
if ¥(X = 1) = p and ¥(X = 0) = 1 - p. 

We say that X has Geom + distribution with parameter a £ (0,1) if 

V(X = k) = a(l - a)*- 1 

for k > 1. If X has Geom + (a) distribution then X — 1 is said to have Geom°(a) distribution. 

Now we define a Bernoulli-geometric distribution, with parameters p and a. A random 
variable with this distribution has the distribution of the product of two independent random 
variables, one with Ber(p) distribution and the other with Geom + (a) distribution. That is, 
A ~ Ber(p)Geom(a) if 

fl-p, fc = 

¥(A = k) = { 

Ipa(l-a)*- 1 , k>\. 

We have Ej4 = p/a, and the probability generating function of A is given by 

•iff ( V A\ _ Cf -p) - (l - P- q> 

1 >~ \-(\-a)z 

The distribution of A, conditioned on being non-zero, is simply Geom + (a). 

In passing, note that such a random variable A may also be represented as a geometric 
number of independent geometries, in the case p < a. Namely, let V ~ Geom°(p), and 
let Wi be i.i.d. Geom + (7), where 7 = (a — p)/(l — p), and independent of V. Then define 
R = Wt + W 2 + ■ ■ ■ + W v . One has R = A. 
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Alternatively, let V ~ Geom (r) where r = (1— p)p/(l — pa), and let Wi be i.i.d. Geom (7), 
where 7 is as above, and independent of V. As before, define R = W\ + ■ • • + Wy 4 , again one 
gets R = A. 



4 Bernoulli-geometric queues 

We consider a queue with arrival process A n , n G Z where A n are i.i.d. Bcr(p)Geom(a), and a 
service process iS n , n € Z independent of the arrival process and with i.i.d. Ber(g)Geom(/3). 
For stability we clearly need ES n > EA n , i.e. p/3 < qa. 

In general, the queue-length process is not reversible. For example, if a is small and (3 
is large, then one tends to see infrequent but large arrival batches, and frequent but small 
departure batches. 

However, under a further condition, which in effect reduces the number of parameters from 
4 to 3, we show reversibility and various other related properties. The relevant condition is 
that 

a P P q 



1-al-p 1 — (3 1 — q 



Note that combined with the stability condition p/3 < qa, condition (4.1 1 implies that (3 < a 
and p < q. 



For a given service distribution, condition (4.1 1 gives one degree of freedom for the arrival 



distribution. In particular, suppose q and (3 are fixed, giving an overall service intensity of 



q/(3. Then for any A < q/(3, there exists a unique pair p, a satisfying (4.1 1 and such that the 
arrival intensity p/a is equal to A. We further discuss the relevance of (4.1 1 to the reversibility 
of the queueing process in Section |4.5| 



Theorem 4.1 Suppose that p/3 < qa, and (4-1) holds 



(i) The queue-length processes (X n ) and (Y n ) are reversible; moreover, they are jointly 
reversible in the sense that 

(. . . , X-i, y_i, x , Yq,Xi,Yi, ...) = (..., x 2 , Yi, Xi, Yq, x , y_i, . . . ). (4.2) 

(ii) The departure process (D nj n S Z) has the same law as the arrival process (A n ,n £ Z). 

(in) The stationary distributions of the queue-length processes (before and after service) are 
given by 

X n ~ Ber(c) Geom (7) , 

where 

(3 1 — a a — (3 

C = \-l3 a ' 7= 1-/3' 

and 

Y n ~ Ber{p + c — pc) Geom (7) . 

(iv) For all n, the queue length X n at time n is independent of the process of departures 
(Di, i < n) before time n. 



(i 



From (iii), one has 



EX 



EY = EX+ P -. 
a(a — p) a 



Proof of Theorem To show the reversibility in (i) and the stationary distribution for X n 
in (iii), it is enough to show that for all fc, to, r, with to > fc and to > r, 

7r(fc)P(r = m\X = k)P(X 1 = r\Y = to) = n{r)F{Y = m\X Q = r)V(X 1 = fc|F = m), 

(4.3) 

where n is the probability mass function for the distribution of X n given in (iii). 

If k = r, this is obvious. The cases fc < r and fe > r are symmetric, and one only needs to 
check one, say k < r. 

Note that 

P(Yo = m\X = k)= ¥(A = m-k), 



while 



F(Xi = r\Y Q = m) = 



F(S>m-r), r = 
F(S = m-r), r>0 



Further one has by assumption that 
ir(x) 



1 — c, x = 

c 7 (l-7) a; - 1 , x>\ 



F{S = x) 



x = 



F(A = x) 



> x) 



I-P, 
pa(l — a) 



x-l 



1, 



x = 
x > 1 

x = 



q/3(l - P) x - X , X>1 \q(l-f3) x -\ r I 

Now one can for example divide into four further cases: (1) k > 0, m = r; (2) fc > 0, to > r; 



(3) fc = 0, m = r; (4) fc = 0, to > r, and check (4.3 1 directly in each case. 

For example, in case (1), we can use the forms of the probability distributions above to 
give 



7r(fc)P(F = m\X Q = k)F(Xi = m\Y = to) 
= ir(k)F(A = m - fc)P(5 = 0) 

= c 7 (l - 7 ) m " 1 (l - p)?/3(l - Z?)" 1 -^ 1 

= 07(1 - 7 ) m - 1 (l - - P)" 1 -"- 1 

= tt(to)P(A = 0)P(5 = m - fc) 

= 7r(TO,)P(r = m\X = m)F(X 1 = fc|F = m) 



p a 1 — q 1 — f3 
1— pi — a q (3 



1-a 

(l-/9)(l-7) 



m — k 



Two lines from the end, we used condition (|4.1| and the fact that ^ ^ = 1 (which 



follows from the definition of 7). 

The other three cases follow similarly and we omit the details. 

The stationary distribution for Y n in (iii) follows from the distribution of X n and the 
fact that Y n = X n + A n with X n and A n independent (for example, multiply the generating 
functions). 
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Given the reversibility, the properties in (ii) and (iv) can be deduced using the same 
arguments that have been used to establish Burke's theorem in various settings (originally by 
Reich for the continuous-time M/M/l queue [TT]). 

First note that D n = Y n — X n+ i and A n — Y n — X n . Hence the reversibility property in 
(i) implies that (_D„)„ e z and (A_ n )„ e z have the same distribution. But A n are i.i.d. random 
variables, so (A_„)„ e z and (A n ) ne % in turn have the same distribution, giving (ii). 

Now note that, by ( |2.3[ ) for example, we can write X n as a function of (Ai)i <n and (Si)j< n . 
Hence X n is independent of (Ai)i> n . But, using (i) again, the collection 

{Xn\ A n: A n+ i, A n+ 2i ■ ■ ■) — {X n ; Y n — X n , Y n+ i — X n+ i, Y n+ 2 — X n+ 2, ■ ■ ■ ) 

has the same distribution as the collection 

(X n ; D n -Xi D n -2, -Dn-3j • • • ) = {X n ', Y n -i — X n , Y n -2 — X n -i, Y n -z — X n+ 2, . . . ) . 

Thus indeed X n is independent of {Di) i<n as required for (iv). □ 



4.1 Fixed points 

Consider a queueing server defined by a given distribution of the service process. We may 
ask for distributions of the arrival process with the following property: when such an arrival 
process is fed into the queue (independently of the service process), the resulting departure 
process has the same law as the arrival process. Such a distribution of the arrival process is 
called a fixed point for the given service process. 

Starting from Burke's theorem for a continuous-time M/M/l queue, questions concerning 
fixed points have been extensively studied in the context of single-server queues such as 
the model of Section 2.1 - see for example [9] and references therein. They have also been 
considered, although less often, in the context of the model of batch arrivals and services in 
discrete time considered here - see for example [3]. 



Part (ii) of Theorem 4.1 is such a fixed point result. Fix a service process of BerGeom 
type, specified by the parameters q and f3. Let /i = q/(3 = E S n be the service intensity. Then 
for each A < /x, there exists an arrival process which is a fixed point of the queue, and which 
has arrival intensity E A n = A (simply choose p and a to satisfy p/a — A along with condition 
(4.1 1 - there is a unique way to do this). 

In fact, Theorem 5 of jl] implies that this gives the unique fixed-point arrival process 
which is ergodic with arrival intensity A. Furthermore, this fixed point is attractive; loosely, 
this means that if any ergodic arrival process is fed into a tandem of queueing servers with 
this service distribution, the distribution of the resulting output process converges to the fixed 
point as the length of the tandem grows. See [1] for precise definitions. 

We note in passing that if (4.1 1 is satisfied, then the distribution Ber(p)Geom(a) has min- 
imal relative entropy with respect to the distribution Ber(q)Geom(/3), out of all distributions 
with mean p/a. See the introduction of [1] for related discussions. 



4.2 Tandems 

Using Theorem |4. 1| we can also describe the behaviour of systems of queues in tandem. 
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Consider a system of R queues in tandem. Each queue has an independent service process, 
Sn ^ for the rth queue, and each of these processes is a collection of i.i.d. Ber(q)Geom(/3) 
random variables. 

The first queue has an arrival process A„ , which is independent of the service processes 
and is a collection of i.i.d. Ber(p)Geom(a) random variables. We assume that p, q, a, (3 



satisfy (4.1l and the stability condition pfi < qa. 



Now, recursively, let the arrival process to the rth queue be given by the departure process 
from the (r— l)st queue, for r = 2, 3, . . . , R; that is, ^ — D„ ■ Thus a customer departing 
from queue r — 1 moves immediately (within the same time-slot) to the next queue. 

Using Theorem |4 . 1 1 and well-known methods, we obtain the following result: 



Theorem 4.2 (i) All the departure processes have the same distribution, which is 

also the distribution of A^ . 

(ii) The vector Xn , X n , • • ■ , X n of queue-lengths of the R queues at a fixed time n is a 
collection of i.i.d. random variables, whose common distribution is that given in Theorem 

Eur**;. 



Proof: Part (i) is obtained by applying Theorem 4.1 ii) repeatedly. The argument to obtain 



the product form result in (ii) from Theorem 4.1 iv) is exactly the same as for the familiar 



case of M/M/l queues in tandem - see for example Section 2.2 of [TJ. □ 
The result also extends easily to cases where the parameters of the service process vary 
between queues: S n ~ Ber(<7 r )Geom(/3 r ). However, all the pairs (q r ,(3 r ) must still belong to 



the same one-parameter family satisfying (4.1) 



(r) 

We could also consider a vector of "queue-lengths before service" , Y„ . Observe that the 
way we have defined the system of queues in tandem, a customer may be present after arrival 
and before service in several different queues at the same time-slot (since a departure from 
queue r — 1 at time n arrives at queue r at the same time n). So in this case, the corresponding 
result is in fact that Yn , Y^_ t , . . . , Y^_ R+1 form an i.i.d. sequence. This may also be proved 
by similar methods. 

4.3 General service-batch distributions 

In order for the reversibility properties in Theorem |4.1[i),(ii),(iv) to hold, one needs the 



condition (4.1 1 relating the parameters of the distributions of arrival and service distributions. 



However, one may wonder whether this is necessary to have a result like Theorem 4.1 'iii) on 
the stationary distribution of the queue-length. 

In fact, such a property holds in a much more general case. We will still assume the same 
sort of arrival process, but now we will consider any service process which has i.i.d. entries 
which are non-negative integers. 

Theorem 4.3 Suppose (S n ) and (A n ) are both i.i.d. sequences taking non-negative integer 
values, and independent of each other, with E S n > EA n . 

(a) If A n has a BerQGeomQ distribution, then both X n and Y n have Ber()Geom() distri- 
butions. 
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(b) If A n takes values and 1 only, then X n has a GeomP distribution. 

(c) If A n has a Geom (respectively Geom + ) distribution, then also Y n has a GeomP (re- 
spectively Geom + ) distribution. 

Parts (b) and (c) are both special cases of part (a) . Part (c) was already observed in Propo- 
sition 12 of |2J, and results like part (b) for M/GI/l queues are certainly well known. 

Proof: We will use a representation of the queue length as the future maximum of a random 
walk. Arguments of this kind are rather classical; for example, see the book of Takacs [T3] 
for many examples (often in the context of queueing theory). 



Using (2.3), we can write 

m m 

X = maxV (S-r + A- r ) , Y = A -hmaxy^ (-S- r + A_ r ) . 

m>0 ^ — 4 m>0 — ' 

r— 1 r— 1 

Consider for example part (b). Xq is the maximum of the walk 

0, + A_ 1: + A_ t - S- 2 + A_ 2 , .... 

The maximum is almost surely finite since the walk has negative drift. Any positive step of 
this walk has size exactly 1 (since the A n are or 1 and the S n are non-negative). Hence any 
new maximum must be precisely 1 greater than the previous maximum. We may treat these 
maxima as renewal times; the future evolution of the walk (relative to its current position) 
has the same distribution as the original walk. Note that Xq > k iff this walk reaches level k 
at some point. Hence we get 

P(A > k + l\X >k)= V(X > 1\X > 0) = P(A > 1) 

for all k. Thus Xq indeed has a geometric distribution as desired (and by stationary the same 
is true of X n for all n). 

For parts (a) we generalise this argument slightly. Now the arrival batches may have size 
greater than 1, but we will regard such a batch as a sequence of individual steps up, each 
of size 1. Since the arrival batches are BerGeom, we have the memoryless property for A n : 
the distribution of A n — k, conditional on A n > k, is the same for all k > 1. We regard 
the walk as a sequence of steps up of size 1, separated perhaps by some number of jumps 
down (which may be 0). By the memoryless property, the jumps down which separate each 
pair of steps up form an i.i.d. sequence, and so we have a renewal property after each step 
up, and thus in particular after each new maximum. Thus P(A"o > k + 1\Xq > k) doesn't 
depend on k, for k > 1, and indeed Xq has a BerGeom distribution. (The difference from 
the previous paragraph is that the distribution for k = may be different; hence we get a 
BerGeom distribution in general rather than the Geom distribution we had in the specific 
case above). The same argument applies also to lo- 
in the special case (c), the arrival batches are Geom + so that Yq must always be strictly 
positive; hence in fact Yq is itself Gcom + . If the arrival batches are Geom , we could add 1 
to every arrival batch and to every service batch to arrive at the previous case; one obtains 
that Yq + 1 is Geom + and hence that Yq has a Geom distribution. □ 
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4.4 Continuous models 

All the results above have equivalent versions in the case where geometric distributions are 
replaced by exponential distributions. 

X is said to have Ber(p)Exp(a) distribution if 



F{X > x) 



1, x = 

pe~ ax , x>Q. 



X may be represented as the product of a Bernoulli random variable and an exponential 
random variable, or as the sum of a geometric number of i.i.d. exponential random variables. 



Now one could write versions of Theorems 4.1 and 4.2 where the Ber()Geom() distributions 



are replaced by Ber()Exp() distributions. These continuous versions can be proved directly, or 
derived from the discrete versions by taking an appropriate limit under which the parameter 



of the geometric distribution tends to 0. The analogous condition in place of (4.1l is that 
ap/(l — p) = (3q/(l — q). Theorem |4.3| extends similarly (and we no longer need to assume 
that the services S n take integer values). 

An alternative extension is to systems in continuous time rather than discrete time. Write 
the Bernoulli parameters as p = eX and q — e/i; now let e — > and rescale time by e. This 
produces a model in which arrival and service batches occur at the times of independent 
Poisson processes, with parameters A and \x respectively. Again, analogous versions of all 
our main results can be straightforwardly formulated. Note that in this case, the distinction 
between queue-length after service and queue-length before service disappears. 

These continuous limits, both in time and in workload, are also exploited in the next 
section in the context of first-passage percolation models. 



4.5 Role of condition (4.1) 



It may be helpful to give some brief comments about the relevance of condition (4.1 1 to the 
reversibility of the queue- length process as in Theorem |4.1| 

Consider a "busy period" of the process, i.e. an excursion away from 0. Suppose that the 
sequence of arrival batches in the busy period is ax, a%, . . . , a n and the sequence of departures 

is d\ , d®, . . . , d n - We have a\ -\ h a n = di H — • + d n , and that a% + ■ ■ ■ + a m > d\ ^ — • + d m 

for m — 1, 2, . . . , n — 1. Note that this also gives in particular that ai > and d n > 0. 

For reversibility of the queue, we need the likelihood of such an excursion to be invariant 
under time-reversal (in which the role of arrivals and departures is exchanged) . This likelihood 
is given by 



^Xo = 0, Ax = oi, . . ., A n = a n ,Di =di,...,D n = d n j, 
which may also be written as 

IP^Ao = 0, Ax = ai, . . . , A n = a„, Sx = dx, ■ ■ ■ , SVt-i = d n -x, S n > d, 

Note that when we translate from the variables Di to the variables Si, the last equality 
becomes an inequality, since it is at this point that the queue-length returns to and so some 
part of the final service S n may be unused. 
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First consider the case where the arrivals are i.i.d. Geom°(a) and the services are i.i.d. 
Geom°( / 9). Then the likelihood above is equal to 

a n (l - a) 01+ - +a "/3"- 1 (l - /3) d i+-+ d ". 

Since the sum of the at is equal to the sum of the di, this is invariant under the exchange 
(ai, . . . ,a n ) <-> (d n , . . . ,d\) as required. 

Now consider the case where Ai ~ Ber(p)Geom(a) and Si ~ Ber((7)Geom(/3). Now the 
likelihood above becomes 



p"a"(l-a) Ql+ - +a " ( ) g "ft"- 1 (l-fl*+-+M 1 g1 / ) 

VP"/ V 9 P J 

Since in general the number of which are zero may be different from the number of di 
which are zero, this likelihood is invariant under the exchange (oj, . . . , a n ) <-> (d n , . . . ,di) 



only if condition (4.1 1 holds. Hence (4.1 1 is also necessary for reversibility of the queue-length 
process. 

By decomposing the queue-length process into its excursions and arguing in this way, we 
could in fact arrive at a proof of Theorem |4 . 1 1 which avoided many of the calculations needed 
in the proof given above (at the expense of complicating the structure of the proof a little). 

The following property is also related to the discussion above. If A ~ Geom°(a), then 
the ratio W(A = k + l)/¥(A — k) is the same for all k > 0, namely (1 — a). If instead 
A ~ Ber(p)Geom(a), then the same is true except fo r k = 0; in the case k = 0, the ratio is 



P 

multiplied by a further factor . Condition (4.1 1 says that this "adjustment factor" 

1 — a 1 — p 

is the same for the distribution of arrivals as it is for the distribution of services. 



5 First-passage percolation models 

In this section we consider various directed first-passage percolation models, and use Theorems 
14. II and 14.21 to calculate exact values for time constants. 

For < (k,l) £ Z 2 , denote by II ((i, j), (k, I)) the set of "directed paths" (i,r\),(i + 
1, rg), ... ,(k, Tk-i+i), where / < r± < T2 < ■ • • < — i + 1 < I. These paths are strictly 
increasing in the first coordinate, and weakly increasing in the second coordinate. See Figure 
[5] for an example. 

Let 5£ be a collection of i.i.d. random variables, and for each path 7 £ H((i,j),(k,l)), 
define the weight of the path 7 by 5(7) = J2( n r)e-y ^n- N° w we define the first-passage time 
from to (k,l) as the minimum weight over all paths from to (k,l): 

F((i,j),(k,l))= rniri 5(7). (5.1) 

T en(0,j), (/=,/)) 

This model has various possible alternative presentations, for example in [ 12] as a first- 
passage percolation model with weights on the edges, where weights on vertical edges are all 
equal to a constant and weights on horizontal edges are i.i.d. 

By Kingman's subadditive ergodic theorem, for any x > 0, there exists a constant f(x) 
such that 

^F((0,0),([xN\,N))^f(x) (5.2) 
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Figure 5.1: An example of a directed path from (1, 1) to (8,6). 



almost surely. 

In [IS], it is shown that if the random variables S r n have Ber(g) distribution, then 



m = <°; , **<'-«>/«■ (5. 3) 

[ WQX - VI - ?) , x>(l-q)/q. 
In [10], related methods are used to show that if the S£ have Geom°(0) distribution, then 



,"• x</3/(l-0), 
l(Vr^3VT+^-l) 2 , x> 0/(1-0), 



while if the have Exp(l) distribution then 



f(x) = (vr+x-iy 



(5.4) 



(5.5) 



The method of [TO] exploits a link between the percolation model and the system of batch 
queues in tandem of the sort described in Section [4~2] Using this method, we can extend the 
results above to the case of a Ber()Geom() distribution. As suggested by the notation, the 
weights in the percolation problem correspond to service times in the queueing model. 

Fix q and and assume that all the are i.i.d. Ber(g)Geom(/?). In this case the service 
rate at each queue is fi = q/0. For any A < /i, we can choose p and a satisfying (4.1 1 and 
such that p/a = A. If the arrival batches are i.i.d. Bcr(p)Geom(a) then the arrival rate is A. 

(Rather than choosing A < /i directly, we may equivalcntly choose a > or choose p < q. 



To express one of the three variables a,p, A in terms of another, we have, as well as (4.1l 
that 



A = - =p 
a 



p 1 — q 1 — 
1-p q 

(1 ~ a)0q 
a 2 (l - - q) + a0q 



1 



Recall that and q are to be regarded as fixed throughout). 

' apply. In particular, the queue lengths Xq , 



Now the results of Theorems 



4.1 



and 



4.2 



are i.i.d., and their expectation is given by 

EXn = 



1 



a 



a 



(5.6) 
(5.7) 



Y R 



(5.8) 
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We may regard this expectation as a function of the arrival rate A and write it as h(X) to 
emphasise this. We also have 



h(X) 



P (1-g) 

(i-p) q 



p(i - g) 

/3(g-p) 



1 



(5.9) 



Now, in Section 3 of |10j . it is shown, by working recursively from representations such as 



( |2.3| ), that 

f> r - sup J ^ A r -F((m,l), (-1,R)) 1 , 

r=l U=m J 

and that by passing to the limit R — > oo and taking a Legendre transform, one obtains 

f(x)= sup {Ax-/i(A)}. (5.10) 

0<A</j 

We apply this result to various cases in turn. 

Example 1. The weights have Ber(q)Geom(/3) distribution. This is essentially the most 
general case of the ones we treat: all the others will be derived by taking some sort of limit 
from this case. 



Using ( |5.10| and plugging in ( |5.6| -( |5.9| , we get 

' p[p(l-q) + (q-p)P] 



f( x ) = 5Z SU P 

P9 pG(0,g) 



(1-P) 



(5.11) 



or alternatively 



f(x) = P sup 



1 



a(l - P - q) + 0q 



P 



(5.12) 



In principle one could solve a quartic equation to put these expressions into closed form, but 
the result is unlikely to be pretty. However, it is straightforward to show that f[x) = if 
x < (1 — q)/q and f(x) > otherwise. 

Example 2. Now consider the case where the common distribution of the weights is 
Ber(g)Exp(l). To obtain this case we can let P — > in the previous example and multi- 
ply by p. We obtain 



/(*) 



sup r 

0<r<l 



qx 



1 



1 — q + rq 1 — r 



Example 3. Now we consider a model where space becomes continous in one direction. 
Consider taking the Bernoulli parameter q to 0, and the space parameter x to infinity. We 
arrive at a model where the first space parameter i is replaced by a continuous parameter. 
For each r, we have a "service process" S r (t). This process is a jump process; events occur 
at times of a Poisson process of rate 1, say, and to each event is associated a weight, which 
is the amount by which the process S r jumps up at the time of the event. These weights are 
independent and each has Geom(/3) distribution. In the queueing model, the weight occurring 
at time t in the process S r corresponds to the amount of service available at queue r at time 
t. 
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The first-passage time can now be defined by 

l 

F((s,j),(t,l))= inf V [S r (u j+1 ) - S r ( Uj )] . 



and the time constants by 



f(y) 



lim 

JV->oo N 



1 P( (0,0), ([Ny\,N] 



(5.13) 



(5.14) 



which, as at (5.2 1, are a.s. constant for each y. 

To obtain this case from ( |5.11 I, we let q — > and set x = y/q. We obtain 



/(y) = /? sup 

aS(/3,l) 



1 



a 



y 



a{\-p) a -(3 



Example 4. We can take the two limits of Example 2 and Example 3 together. Now the 
service processes will consist of weights occurring at times of a Poisson process, with each 
weight having Exp(l) distribution. In this case we get 

1 



sup r 

0<r<l 



V 



1 



In this final case the translation into closed form is more reasonable; one only needs to solve 
a quadratic equation, to obtain 



f(y) 




8y 3 + 20y 2 + 7y + 1 - ^/8y~+l (8y 2 + 3y + l) 



y < i 
y> l. 



Taking certain other limits or special cases in these examples recovers previously known 
results. Taking q = 1 — (3 in (5.11) gives the case of Geom(/3) weights and leads to ( 5.4 1 , 
and further taking (3 — > and multiplying by (3 gives the case of Exp(l) weights and leads 
to (5.5). On the other hand, taking (3 — ► 1 in (5.111 gives Bernoulli(g) weights and leads to 

(H. 

One could also consider the continuous model of Examples 3 and 4 in the case where each 
weight has value 1, so that the service processes are simply Poisson processes. To do this we 
take /3 — * 1 in Example 3, to obtain simply 

/(y) = ([^-i] + ) 2 - 

Finally, by taking an appropriate continuum limit in any of these cases, one can arrive 
at the Brownian first-passage percolation model which has been quite widely studied (see for 
example pQ, [5]). 
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